Electronic phase separation in the pyrochlore double-exchange model 
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Electronic phase separation and related inhomogeneity is ubiquitously seen in strongly-correlated 
systems. A typical example is found between ferromagnetic metal and antiferromagnetic insulator 
in CMR manganese oxides. Here we demonstrate that the geometrical frustration brings distinctive 
aspects into the phase separation phenomena. From Monte Carlo simulation and a simple energy 
comparison for the pyrochlore double-exchange model, we show that such phase separation takes 
place between ferromagnetic and paramagnetic metals. We discuss the relevance of our results to a 
spin-glassy metallic phase found in Mo pyrochlore oxides under external pressure. 

PACS numbers: 71.30.+h, 71.27. +a, 71.20.Be, 71.10.Fd 



Electronic inhomogeneity is one of the key concepts 
in strongly correlated electron systems. It emerges from 
keen competition between different electronic phases, and 
gives rise to fascinating phenomena, such as nanometer- 
scale domain structure and gigantic response to exter- 
nal perturbations. Many intriguing topics related to the 
electronic inhomogeneity have been studied intensively, 
such as granular superconductivity in cupratesll[ , colos- 
sal magneto- resistance (CMR) in manganites [2[ , and re- 
laxor response in dielectric compounds |3(. 

One of the most interesting, well-known examples is 
found in CMR manganites. It was argued that the com- 
petition between ferromagnetic metal (FM) and antifer- 
romagnetic insulator (AFI) results in an electronically- 
inhomogeneous state in the presence of quenched disor- 
der, and that the inhomogeneous state is highly sensitive 
to external magnetic field and gives rise to the CMR ef- 
fect 0| . Theoretically, the double-exchange (DE) model 
and its extensions have been extensively studied, and it 
was clarified that the model has an instability toward the 
electronic phase separation (PS) The PS occurs with 
a discontinuous change in the electron density, e.g., be- 
tween FM and AFI: The former appears in a wide range 
of density by the DE ferromagnetic interaction of kinetic 
origin 0] and the latter is stabilized at a commensurate 
electron density by the AF super-exchange (SE) interac- 
tion [(| The energetically-competing phases lead to 
the electronic PS as well as bicritical behaviors, which 
are discussed as key players in the fascinating properties 
of CMR manganites. 

In this Letter, we elucidate a new type of electronic PS 
in the DE model, which emerges when the lattice struc- 
ture is geometrically frustrated. Here we show that PS 
takes place between FM and a paramagnetic state for the 
model defined on the pyrochlore lattice which has a three- 
dimensional network of corner-sharing tetrahedra [inset 
of Fig. QJa)]. The paramagnetic state does not have a 
charge gap, in contrast to the AFI state in the unfrus- 
trated cases. While the PS was observed in the previous 
study Q, the present work elucidates its ubiquitousness 
and origin. It should be noted that this PS with a dis- 



continuous change of the density is non-trivial, since the 
transition from ferromagnet to paramagnet is continuous 
in general. This is in sharp contrast to the unfrustrated 
case where magnetic transition between ferromagnet and 
antiferromagnet is necessarily of first order, which can 
naturally accompany a jump in the electron density. 

In addition to the new theoretical aspect, there is an 
experimental motivation for the frustrated DE models. It 
was proposed that the DE mechanism is relevant to FM 
state observed in Mo pyrochlore oxides i? 2 Mo 2 07 (R is 
rare earth) [9]. Recently, it was found that the com- 
pounds exhibit puzzling behaviors under external pres- 
sure [lOj : FM is suppressed and taken over by a peculiar 
diffusive paramagnetic metal (PM). In the intermediate 
region, the system shows a spin-glassy metallic (SGM) 
behavior, which is in contrast with a spin-glassy insulat- 
ing behavior observed in CMR manganites. The origin 
of SGM is not clear yet. We will discuss the relevance of 
our results to this SGM state. 

Our model is the DE model defined on the pyrochlore 
lattice, shown in the inset of Fig. []Ta) . We consider the 
transfer integral and the SE coupling for nearest-neighbor 
(n.n.) sites only (the next n.n. hopping will be included 
later). For simplicity, we take the limit of large Hund's- 
rule coupling and treat the localized spins Sj as classical 
vectors with a normalized length Si = 1. Then the 
model is effectively described by a spinless Hamiltonian 
in the form [H El 

u = E [ - '« & cj + h.c.) + J AF Si ■ Sj] , (1) 

where the first term describes the hopping of electrons 
whose spins are locally forced to align with the localized 
spins, and the second term represents the AF SE inter- 
action between the localized spins. Here the effective 
transfer integral is given by Uj — t [cos(0,/2) cos(9j/2) + 
sin(9i/2)sm.(9j/2)exp{—i((j)i — <f>j)}], where (6i,4>i) de- 
scribes the angle of the localized spin S^. We set an en- 
ergy unit as t = 1 and the Boltzmann constant fcs = 1. 

We calculate thermodynamic properties of the model 
([T]) by Monte Carlo (MC) simulation (gj, for the system 
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FIG. 1: (Color online) (a) Monte Carlo results for the electron 
density n as a function of the chemical potential /i at Jaf = 
0.04 and T = 0.012. The inset shows the cubic unit cell of the 
pyrochlore lattice, (b) Spin correlation (Si • Sj) for two states 
on the verge of the phase separation at /Lt c i — 1.4, plotted as a 
function of the distance in unit of the nearest-neighbor bond 
length. The data are measured along the ID chains in the 
pyrochlore lattice, (c) and (d) Density of states for the two 
states. The data in (b)-(d) are for iV s = 128 sites and the 
errorbar is in the symbol size or the line width. 



sizes N s = 16, 32, 64, and 128 sites which correspond 
to lxlxl, 1x1x2, 2x2x1, and 2x2x2 cubic 
unit cells of the pyrochlore lattice shown in the inset of 
Fig.QJa). To suppress the finite-size effect, we take an av- 
erage over the twisted boundary conditions p^| - [l5| which 
is a well-established technique to smear out the extrin- 
sic "staircase-like" structure of the electron density as a 
function of the chemical potential. This averaging is effi- 
cient, in particular, to identify PS region, since it enables 
us to extract the intrinsic jump associated with PS even 
in small-size clusters, as demonstrated below. 

Figure QJ a) shows a typical MC result for the electron 
density per site n = ^i{c\ci) /N s as a function of the 
chemical potential /j,. It is clearly seen that the system- 
size dependence is suppressed even at the low tempera- 
ture (T) , owing to the averaging over the twisted bound- 
ary conditions. The results reveal two sharp jumps: One 
is from n 0.49 to n ~ 0.75 at fx c \ ~ 1.4, and the other 
is from n ~ 0.04 to n ~ 0.16 at fi C 2 ~ —3.2. These two 
jumps signal the electronic PS; the system is not stable 
at a density in these regions and phase-separated into 
two states at the edges of jump. 

Both the PS take place between ferromagnetic and 
paramagnetic states. This is corroborated by the spin 
correlation plotted in Fig. QJb) . The results show typical 
behaviors on two sides of the PS at /i c i ~ 1.4. At fi < fi c i, 
the spin correlation is large and converges to a positive 
value for farther neighbors, indicating a long-range fer- 
romagnetic ordering. On the other hand, at \i > fi c i, the 
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FIG. 2: (Color online) Phase diagrams determined by MC 
calculations at (a) Jaf = 0.00, (b) 0.04, and (c) 0.08. Symbols 
surrounding the PS regions show the critical densities where 
n changes discontinuously as a function of /i. Filled circles 
at T — correspond to the phase boundaries in Fig. EJb). 
Crossed squares indicate the second-order phase boundary be- 
tween FM and PM. The lines are guides for the eye. 



spin correlation is AF for neighboring sites, but decays 
quickly to zero, suggesting a disordered paramagnetism. 

The ferromagnetic and paramagnetic states are both 
metallic. Figures [IJc) and []Jd) present the density of 
states (DOS) per site for the two states. The results 
show that both states have no energy gap at the Fermi 
energy uj = 0. This is consistent with the fi dependence 
of n in Fig. [Ha): dnjd[i is nonzero on the both sides 
of PS, namely, the system is charge compressible. The 
ferromagnetic state is stabilized by the DE interaction, 
and hence, it is metallic. The paramagnetic state is also 
considered to remain metallic; a gap formation by a sim- 
ple AF ordering is suppressed by strong frustration [l6| . 
All the above situations of the magnetic and electronic 
states hold in a similar manner for PS at a lower density 
at [i c i ~ —3.2, while in this case the system is the ferro- 
magnetic (paramagnetic) state for fj, > /i C 2 (fx < ^02)- 

We therefore conclude that the electronic PS occurs 
between the FM and PM states. This is a new type of 
PS appearing between metallic states, in sharp contrast 
with that in the unfrustrated case between FM and AFI 
with a clear energy gap [g, [7[ ■ In the present model, we 
have shown in our previous paper Q a crossover from 
the simple PM to the cooperative PM (CPM) states; the 
crossover temperature generally goes to zero on the verge 
of PS because of the cancellation between DE ferro and 
SE AF interactions. Therefore, the electronic PS dis- 
cussed here should in general occur above the crossover 
temperature and thus between FM and the simple PM. 

The finite-T phase diagram for PS is determined by 
keeping track of the discontinuity of n(/u). The results 
are shown for the higher-density PS at /i cl for different 
values of Jaf in Fig.O The PS region appears in a dome- 
like shape in the plane of n and T, whose top connects to 
the second-order phase boundary between FM and PM. 
As increasing Jaf, the maximum temperature of PS is 
first enhanced slightly, but is suppressed for larger Jaf: 
PS disappears for Jaf > 0.1 in the T range that we have 
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FIG. 3: (Color online) (a) Comparison of the ground- 
state energies between the completely-disordered PM and the 
perfectly-ordered FM. The dashed lines connecting the circles 
are common tangent lines for two energy curves, identifying 
the PS regions. The insets are the density of states for PM 
(upper inset) and FM (lower inset), (b) Phase diagram de- 
termined by the energy comparison. See the text for details. 



calculated in MC (T > 0.01). 

Now we try to understand the origin of PS on the ba- 
sis of a simple argument on the energy comparison. Here 
we assume that FM is in the perfectly-ordered (fully- 
saturated) state and PM is in the completely-disordered 
state which is approximately obtained by using PM solu- 
tion of the dynamical mean- field theory [17j . In this as- 
sumption, we neglect not only the short-range AF corre- 
lations in the PM states, but also possibility of any other 
symmetry-broken states. We will see, however, that the 
energy comparison between the two simple states gives 
a reasonable estimate of PS regions in comparison with 
the fiuite-T MC results. 

Figure [21a) shows the results for this analysis. The 
lowest curve in the figure is the energy for the FM state 
at Jaf = 0, which is essentially given by that for non- 
interacting spinless fermions. The FM state has two dis- 
persive bands for —6 < u> < — 2 and —2 < u < 2, and 
two degenerate flat bands at u = 2. Because of the 8- 
functional flat bands, the energy as a function of n is 
given by a straight line from n = 0.5 to n = 1.0 in 
Fig. EJa). DOS is shown in the inset with the one for the 
PM state. When we turn on Jaf, the perfectly-ordered 
FM state costs the SE AF energy, though the completely- 
disordered PM state does not because the SE energy is 
averaged to be zero. Figure^ a) illustrates the situation 
for Jaf = 0.04 and 0.08. We can identify the PS region 
by drawing common tangent lines for the energy curves 
of FM and PM states (dashed lines in the figure). Note 
that PS appears in both small and large density regions 
for Jaf > 0. 

We carry out the energy comparison with varying Jaf, 
and summarize the phase diagram shown in Fig. [3Jb) . In 
the middle of the phase diagram, we have the DE FM 
state, which shrinks as increasing Jaf because of the loss 
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FIG. 4: (Color online) Phase diagram determined by the 
energy comparison for the model with next nearest-neighbor 
hopping t' at (a) t' = 0.1 and (b) 0.2. 



of kinetic energy. The FM state is surrounded by the PS 
region between FM and PM [l]|. The PS has a finite 
width in all regions and a direct continuous transition 
from FM to PM is prohibited, except for a marginal point 
at (n, J AF ) ^ (0.3,0.11). 

Let us compare the energy argument and the finitc- 
T MC results. The phase boundaries obtained from the 
analyses in Fig. [3] are plotted as the filled circles at T = 
in Fig. [2] The finite-T MC phase boundaries appear to 
be extrapolated smoothly to the T = results, except at 
Jaf = 0, where some discrepancy in the right boundary is 
seen. It is probably due to the singular ^-functional DOS 
discussed in Fig. (3Ja) . The overall agreement supports 
that our simple energy comparison is not far from the 
reality and that the PS is essentially a metal-to-metal 
one from FM to PM. 

This peculiar PS is not specific to the present model 
(fT]). We demonstrate it here by extending the model with 
including the next n.n. hopping t' . Figure 3] presents the 
results of the similar energy comparison for the extended 
model. For finite t' , the flat-band singularity in FM is 
removed, at the same time, the sharp cutoff at the band 
top in PM is rounded off [cf. the insets of Fig. |3Ja)]. 
Nonetheless the metal-to-metal PS remains robust in a 
similar manner, as shown in the phase diagrams in Fig. [4] 
PS will be robustly observed even when the model is 
extended to describe more complicated band structure 
and orbital degeneracy. 

We compare our results with those obtained in the ab- 
sence of frustration. In the absence of frustration, PS 
takes place between FM and AFI as mentioned in the in- 
troduction. AFI appears at a commensurate filling, such 
as n — 1.0. The bandwidth of the AFI state is much 
smaller than that of FM, which drives the discontinuous 
change of the electron density associated with PS. Fur- 
thermore, the discontinuity is naturally expected since 
the magnetic transition between ferromagnet and anti- 
ferromagnet is necessarily of first order. In contrast to 
the unfrustrated case, our PS obtained for the frustrated 
pyrochlore model occurs between FM and PM. The se- 
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vere frustration suppresses the AF long-range order. In 
the ground state, the DE-induced ferromagnetism bears 
full saturation of moment to maximize the kinetic energy, 
and hence, if there is a transition from FM to a param- 
agnet by changing the chemical potential, the transition 
becomes discontinuous with a jump of the moment. The 
magnetic discontinuity naturally accompanies a discon- 
tinuous change of the electron density, which leads to 
PS, as demonstrated in Figs. H and H The AF SE in- 
teraction enlarges the paramagnetism to doped metallic 
regions with reducing the FM region, resulting in the pe- 
culiar metal-to-metal PS in our model. Hence the com- 
plete polarization of the DE FM state at T = plays a 
key role in the present PS. 

We note that the argument is in some sense precipitate 
and the situation will be more complicated in the low-T 
limit for large Jaf region, since the spin disordered state 
might be unstable toward some symmetry breaking, such 
as an incommensurate magnetic ordering. It is antici- 
pated in the present DE model that the kinetic motion of 
electrons, which in general leads to farther-neighbor DE 
interactions, acts as a degeneracy-lifting perturbation. In 
fact, in the pyrochlore Heisenberg AF spin system, it is 
known that any perturbation such as farther-neighbor 
exchange interactions may force the system to order . 
What we have revealed is that the T scale for such an- 
ticipated phase transitions is very small in the present 
DE model (at least, smaller than the T range reached 
in the MC simulation) and the PS between FM and PM 
takes place well above it. We expect that a similar PS is 
widely seen in the DE systems on other frustrated lattice 
structures. 

Let us discuss the relevance of our results to experi- 
ments. Mo pyrochlore oxides -R2M02O7 exhibit peculiar 
phase transitions under external pressure [l(|. FM be- 
comes unstable and is taken over by a diffusive PM; in be- 
tween, a spin-glassy metallic (SGM) state appears at low 
T. In the transition to SGM, there is a glassy response in 
the AC susceptibility but no substantial anomaly in the 
resistivity. This SGM is peculiar because a spin-glassy 
behavior is observed in the insulating phase in the case 
of CMR manganites. The spin-glassy insulating behav- 
ior can be understood by competition between FM and 
AFI with a formation of FM clusters in the matrix of 
AFI (or charge/orbital ordered insulator). In the Mo py- 
rochlores, the competing phases under pressure are FM 



and PM, and therefore, we speculate that the interven- 
ing SGM state originate from the electronic PS found in 
our results. The glassy magnetic response is presumably 
explained by domain-like structure or a mixed state; for 
such situation, it is crucial to consider the effect of long- 
ranged Coulomb interaction as well as quenched random- 
ness [13,[2l|. In the present case, the competing phases 
are both metallic, and hence, we expect that the resulting 
glassy state remains metallic. Further study by extend- 
ing our model is necessary to confirm this scenario. It is 
also interesting to extend the experimental study away 
from the commensurate filling for comprehensive under- 
standing of the phenomena. 

Finally we comment on the relation between PS and 
the flat band singularity at Jaf = 0. It was argued 
that the DE FM is no longer stable when the chemi- 
cal potential is in the flat band (22|. This instability 
corresponds to the jump from FM at n = 0.5 to the 
macroscopically-degenerate state at n — 1.0 in Fig. |3l A 
jump of electron density associated with the flat band 
singularity smoothly evolves to PS for Jaf > 0. This 
suggests a possibility to understand the instability from 
the viewpoint of the electronic PS. Further analysis is left 
for future study. 

To summarize, we have investigated the electronic 
phase separation in the double-exchange model defined 
on the frustrated pyrochlore lattice. The phase separa- 
tion is found to occur between the ferromagnetic metal 
and the paramagnetic metal: The former is stabilized by 
the double-exchange interaction, and the latter is induced 
by the frustrated antiferromagnetic super-exchange in- 
teraction. We discuss that the phase separation ex- 
plains the peculiar spin-glassy metallic behavior in Mo 
pyrochlore oxides under external pressure. The metal-to- 
metal phase separation is characteristic of the frustrated 
double-exchange system, which is not seen in the unfrus- 
trated models studied for the colossal magneto-resistive 
manganites. 
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